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SUMMARY 


I' 

For  a  nuR'.ber  of  turbulent  flovj  problems,  a  first-order  closure 
or  eddy  transport  approach  to  the  modellnn  of  the  turbulence  is 
not  sufficiently  accurate.  It  is  necessary  that  a  more  powerful 
approach  be  used  to  study  the  dynamics  of  the  turbulence.  Flows 
involving  chemical  reactions,  such  as  those  in  a  chemical  laser, 
are  good  examples  of  flows  requiring  a  more  powerful  approach.  The 
mixing  in  HF  and  DF  chemical  lasers  takes  place  under  conditions  of 
very  large  heat  release.  The  coupling  between  the  heat  release  and 
the  turbulence  is  an  important  feature  of  the  flow.  A.R.A.P.  has 
been  funde^tm'd'er' APOSR  Contracts  F4il620-73-f^-0027  and  P446-20-75-C- 
0926^to  employ  this  basic  method  to  investigate  problems  connected 
with  chemical  lasers  of  particular  Interest  to  the  Air  Force  Weapons 
Laboratory.  These  studies  have  resulted  in  considerable  insight 
into  the  complex  processes  occurring  in  laser  cavities.'^  These 
studies  have  been  documented  in  References  1,  2,  and  3- 

■'This  report  describes  recent  results  of  .studies  using 

the  reacting  shear  layer  (RSL)  computer  program  and  discusses  the 
proposed  procedure  for  handling  multi-step  chemical  reactions  in  a 
second-order  closure  computation  of  turbulent  reacting  flews. ^ 
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I .  TNTRODUCTIOn 


A  second-order  closure  program  for  the  lnvestir,ation  of  tur¬ 
bulent ,  compressible,  chemically  reacting  flowfields  is  under 
development  at  A.R.A.P.  Second-order  closure  methods  have  provided 
an  understanding  of  the  physics  of  tui’bulent  flov/s  which  could  not 
be  provided  by  classical  first-order  closure  techniques  such  as 
the  eddy  viscosity  method.  We  have  successfully  applied  these 
techniques  to  a  variety  of  turbulent  flov/  problems  includinr:  iricom- 

pressible  and  compressible  shear  layers, the  behavior  of  the 

4 

planetary  boundary  layers,  and  the  dist^eiusal  of  chemically  reacting 

q 

pollutants  in  a  turbulent  atmosphere.' 

A.R.A.P,  is  funded  under  AFOSR  Contract  Ho.  li6<?0-7  j-C-0026 
to  develop  an  invariant  second-order  closure  model  capable  of  investi¬ 
gating  the  mixing  and  reactions  in  HP  and  DF  chemical  lasers.  The 
mixing  in  these  devices  takes  place  under  conditions  of  large  Iccal 
heat  release  due  to  the  highly  exothermic  chemical  reaction.  The 
interaction  between  the  heat  release  and  the  turbulence  is  an  impor¬ 
tant  feature  of  the  flow. 

The  set  of  second-order  closure  equations  governing  the  mixing 
and  reaction  of  compressible,  turbulent,  multi-component  reacting 
shear  layers  has  been  modeled  and  programmed  for  solution  in  our 
RSL  program.  This  interim  scientific  report  describes  some  recent 
results  of  studies  using  this  prorram  and  discusses  the  proposed 
procedure  for  extending  the  program  for  studying  multi-step  chemical 
reactions  in  a  second-order  closure  computation. 


reacttih:  shear  layer  (RSL)  PROHRAH 


II . 

The  reacting  shear  layer  program  has  been  completely  checked 

out  ever  the  past  six  months.  Program  predictions  for  nonreacting 

6  T 

flows  ’  (two-dimensional  wake,  heated  wake,  heated  planar  jet) 
are  in  good  agreement  with  experimental  measurements.  Comparison 
of  program  predictions  with  experim.ental  data  for  reacting  flovjs 
is  now  underway.  The  RSL  program  has  been  used  for  computations 
of  the  HP  chemical  laser  flovjfield,  hydrogen-air  premixed  combustion 
and  diffusion  flames  and  propane-air  premdxed  com.bustion  and  dif¬ 
fusion  flames.  The  calculations  have  been  performed  using  tv;o  dif¬ 
ferent  models  for  the  scalar  correlations.  One  model  (used  during 
the  model  development  stage)  was  to  set  all  third-order  (and  higher- 
order)  scalar  correlations  to  zero.  The  other  model  is  the  A.R.A.P. 
"typical  eddy"  box  model.  Preliminary  results  for  exothermic  react- 
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ing  flows  suggest  that  there  is  only  a  small  difference  between  the 
calculations  using  these  two  m.odels.  More  extensive  tests  are 
necessary  to  check  the  validity  of  this  result  for  different  flow 
geometries  and  in  different  flowfield  regions.  It  is,  of  course, 
obvious  that  if  the  model  involving  setting  the  higher-order 
correlations  to  zero  can  be  used  for  most  of  the  calculations,  it 
v;lll  lead  to  considerable  program  simplification  and  increase  in 
speed . 

The  use  of  the  "typical  eddy"  box  model  in  these  calculations 
has  verified  the  basic  concept  of  modeling  the  scalar  correlations 
using  a  probability  density  function  composed  of  a  set  of  Dirac 
functions  of  variable  strengths  and  position  in  the  scalar  phase 
space.  The  computations  have  also  allowed  us  to  determine  some 
shortcomings  of  the  present  model.  A  number  of  sample  modifications 
have  already  been  made  to  the  model  to  extend  its  applicability  to 
more  severe  environments ,  and  some  othei-  changes  are  currently  under 
study . 


The  curretit  version  of  the  prof^rarri  Ir  restrict-ed  to  one-step 
forward  reactions.  In  tne  past  month,  we  have  developed  a  procedure 
for  handlinn  multi-step  reaction  systems  inGludint:  the  complete 
turbulence-chemistry  interaction.  The  procedure  l.s  briefly  described 
in  the  next  section. 
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III.  SF.COND-ORDER  CLOSURE  FOR  MULTI-STEP 
REACTION  SYSTEMS 

The  "typical  eddy"  box  model  has  been  most  extensively  used 
up  to  now  with  tf^ree  species  (a,  3  and  y)  flov;fields.  In  principle, 
the  model  is  quite  r.eneral  and  can  be  extended  to  a  larpe  number 
of  species.  However,  the  direct  addition  of  each  new  species  in 
the  second-order  closure  system  will  add  ten  eo.uations  for  the 
means  and  second-order  correlations  and  does  not  appear  to  be  a 
reasonable  procedure  for  complex  reaction  systems.  Actually,  such 
a  procedure  is  not  necessary. 

Consider  the  mixing  and  reaction  of  two  initially  separate 
sti'eams  o  and  B  .  These  streams  can  be  known  mixtures  of  a 
large  numiber  of  chemical  species.  The  streams  mix  and  react 
following  a  large  number  of  elementary  reactions  and  form  a  mixture 
of  product  species  designated  as  y  .  We  assume  that: 

1.  Th.e  composition  of  mixtures  a  and  B  is  fixed  for 
the  entire  flowfield.  There  are  no  internal 
reactions  among  the  species  in  these  two  reactant 
mixtures . 

2.  All  the  product  species  that  compose  y  are 
molecularly  mixed.  There  are  internal  reactions 
within  the  y  species  mixture.  Due  to  the  dif¬ 
ferent  rates  of  tlie  various  elementary  reactions 
(between  species  present  in  a  and  B  mixtures 
and  the  species  present  in  y),  the  composition  of 
the  mixture  y  varies  at  different  points  across 
the  flowfield. 

Consider  the  following  two-reaction  system  as  an  illustration 
of  the  procedure.  The  complete  set  of  reactions  that  we  will  be 
using  in  our  study  of  the  DF  chemical  laser  is  presented  in  Table  I. 
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a  known 

B  iie  +  F  +  F.,  knovm  proportion;’ 

t— 

y  kie  +  +  F.,  +  D  +  F  +  I^F 


The  interaction  between  turbulence  and  chemistry  only  has 
to  be  taken,  into  account  for  reactions  between  a  and  B  , 
ot  and  y  ,  and  6  and  y  .  Species  within  y  are  assumed  to 
be  molecularly  mixed  and,  therefore,  sj'ecies  correlations  do  not 
have  to  be  considered  for  the  internal  reactions.  For  conven¬ 
ience,  we  rewrite  the  reaction  system  ns  piven  below,  tappinp 
each  cliemical  species  wltli  the  mixture  that  it  is  a  part  of: 

D,,(a)  +  F(B)  DF(y)  +  D(y) 

D.,(a)  +  P(y)  DF(y)  +  P(y) 

L 

D^iy)  +  F(P)  ^  PF(y)  -t  D(y) 

D^Cy)  F(y)  -*■  DFCy)  -t  n(Y)  internal  reaction 
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!'’,(rO  +  1'(y)  +  !•’(>) 

F,(y)  +  IHy)  PFCy)  +  F{y)  i iitei’tial  ro'ictif'n 


\\o  oan  now  write  t’ne  cheir.i.'al  oonree  tei’ma  )’oi’  tfie  component 
sreclep.  For  example, 


dP,,(rt) 

dt 


=  -  k- 


-  k. 


P.,(a)F(B)  +  p;(rt)^''(B) 


P.,(a)F (y)  +  ’  (y) 


c1L^(y) 

dt 


P.(y)F(B)  +  p.’ {y)F’  (1^ 
fk;(Y)F(Y) 


,  et  c  . 


Two  other  pieces  of  informatioti  are  now  required  to  co!:ir]ote 
this  foritula t ion .  First,  we  will  he  selv.inr  equations  fe:’  the 
second-order  eerrela  t;  i  ons  a'B'  ,  a'Y'  >  'Tn'd  P'y'  .  ?Jow, 
pT(a)F'(B)  ,  Pi(a)F'(Y),  etc.,  can  be  simply  calculated  from 

c. 

these  correlations  when  one  knows  the  cc^mpos  1 1 i oii  of  the  mixtures, 
and  so  the  source  terms  are  completely  detoimiineii .  Peconci,  we  have 
to  find  the  total  loss  of  a  and  P  .  In  ttie  above  examnlo,  P 
is  beinr  lost  due  to  the  consumption  of  both  F(R')  and  F,^(P'). 

We  determine  tlie  term  tb.at  causes  the  prcater  loss  of  B  and  use 
the  figure.  If  at  some  point  the  consumption  of  F(B)  leads  to 
the  larper  loss  of  P  ,  then  one  lias  to  keep  track  of  the  corres¬ 
ponding  components  F^(P)  and  ne(B)  ttiat  now  become  part:  of  y  . 
The  composition  of  Y  has  to  be  calculated  nl  enc!i  point  across 
the  flowfield. 

This  procedure  operates  within  the  frameworl-:  of  Miree  over¬ 
all  species  mixtures  «  ,  P  ,  and  Y  and  the  models  and  equations 
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,ioV'’ I t.‘ i  e-ii'!!Oi’  can  c-till  be  in'od.  onl;;  chanfa; 

I : ■  i r’.  t !  1  >;■  r. t''  1  u t  :  n  of  t li e  in e a n  c b. cn i c a.  1  r- p e c  i e s  v,' h e r’ ( >  1 
complete  act  of  elementary  reactions.  .At 
rrc'scn*:  t  i  -le ,  v.’e  are  maVrinr,  t)!e  prorr’am  chani'os  to  otn 
c!ien:!cal  ]  aoer  flcwfleld  urinr.  al]  the  clonn'ntary  i 
ILo.t.t'J  in  Table  I. 
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